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The largest Lyapunov exponent of a system composed by a heavy impurity embedded in a chain 
of anharmonic nearest-neighbor Fermi-Pasta-Ulam oscillators is numerically computed for various 
values of the impurity mass M. A crossover between weak and strong chaos is obtained at the same 
value e,, of the energy density e (energy per degree of freedom) for all the considered values of the 
impurity mass M. The threshold coincides with the value of the energy density e at which a 
change of scaling of the relaxation time of the momentum autocorrelation function of the impurity 
ocurrs and that was obtained in a previous work [M. Romero-Bastida and E. Braun, Phys. Rev. E 
65, 036228 (2002)]. The complete Lyapunov spectrum does not depend significantly on the impurity 
mass A/. These results suggest that the impurity does not contribute significantly to the dynamical 
instability (chaos) of the chain and can be considered as a probe for the dynamics of the system to 
which the impurity is coupled. Finally, it is shown that the Kolmogorov-Sinai entropy of the chain 
has a crossover from weak to strong chaos at the same value of the energy density that the crossover 
value Ej, of largest Lyapunov exponent. Implications of this result are discussed. 
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I. INTRODUCTION 

During the past two decades or so, there has been grow- 
ing evidence of the connection between the underlying 
chaotic microscopic dynamics of many-particle systems 
and its observed macroscopic behavior. For example, the 
largest Lyapunov exponent (LLE) Ai, which measures 
the exponential rate of divergence of two originally close 
trayectories in phase space, has been found to be an in- 
dicator (order parameter) of phase transitions Q, |^ . 
Moreover, the largest Lyapunov exponent has also been 
related to the transport coefficients of fluid systems with 
continuous potentials For some special cases, 

e.g., hard-sphere systems, the theory of Lyapunov expo- 
nents is remarkably developed 0- However, considering 
randomness and transport in general as signatures of mi- 
croscopic chaos raises subtle and fundamental issues in 
statistical physics that have to be carefully discussed and 
investigated. For the case of Brownian motion the anal- 
ysis of simplified models such as a labeled (test) particle 
immersed in a one-dimensional system of hard rods 1^ 
and an impurity in a harmonic crystal has shown 
that the motion of the tracer particle may be Brownian 
even when the full dynamical system, particle plus fluid, 
is not chaotic. These examples show that microscopic 
chaos is sufficient, but not necessary to produce Brown- 
ian motion. The main problem then is to prove beyond 
doubt if it is possible to detect any particular feature of 
the microscopic dynamics, either regular or chaotic, in 
the behavior of the tracer particle. Only then it would 
be possible to assess the relevance of the microscopic dy- 
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namics on macroscopic behavior. 

The FPU model, which is a one-dimensional chain of 
nearest-neighbor anharmonic oscillators, is a system that 
has been extensively studied over the past decades with 
relation to the problem of energy equipartition. From 
the dynamical system's perspective, it was the starting 
point in the study of chaotic dynamics in many-degrees- 
of- freedom systems (for a recent review, see Ref. llOj). 
In particular, numerical simulations jllj revealed a rich 
phase space dynamics that is controlled by the energy per 
degree of freedom e = E/N. Two qualitatively different 
regimes exist in the dynamical behavior of the system: 
it is strongly chaotic and phase-space diffusion (to be re- 
ferred from now on as microscopic diffusion) is fast when 
the energy density e exceeds a threshold e^, whereas it 
is only weakly chaotic (i.e., almost periodic) and micro- 
scopic diffusion is slowed down when the energy density is 
below the threshold , which in Ref. Illj] was called the 
strong stochasticity threshold (SST). The detection of 
this transition between two different dynamical regimes 
is performed by means of the LLE Ai(e), which exhibits 
a change in its scaling behavior precisely at the value 
of the SST. 

Recently it was shown that, if the FPU chain is cou- 
pled to a heavy impurity, the latter performs Brownian 
motion [T^ . This system, being one-dimensional, can 
be much easily studied than three-dimensional systems 
with continuous potentials, and offers a convenient start- 
ing point for a systematic study of the relationship be- 
tween the microscopic dynamics and the macroscopic, 
statistical behavior. In particular, it is a suitable model 
to explore the possibility to find some indication of the 
known microscopic dynamics in the statistical behavior 
of the heavy impurity. The results of the present work 
suggest that this is the case indeed. In Sec. II we de- 



scribe the model to be used. In Sees. Ill and IV we 
report the resuhs of the statistical and dynamical behav- 
ior of the system, respectively. Section V is devoted to 
discuss the relationship between the results found in the 
two previous sections. In Sec. VI some conclusions are 
drawn. 



II. THE MODEL AND ITS NUMERICAL 
INVESTIGATION 

In terms of dimensionless variables, the Hamiltonian 
of the model considered in this work is 
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with TOi = 1 for i ^ and toq — M; periodic boundary 
conditions are assumed (x(jv/2)+i = X-m/2)- The model 
describes a system of one-dimensional N coupled non- 
linear oscillators of unit mass with nearest-neighbor in- 
teractions, displacements {xi\, momenta {pi\ and a cen- 
tral oscillator (impurity) of mass M with displacement 
xq = X and momentum po = P. The value P — 0.1 
was used in the computation of all the numerical results 
hereafter reported. This model will be referred to as the 
Modified FPU (MFPU) model. 

As initial conditions we choose the equilibrium value 
of the oscillators displacements, i.e. Xi{Q) = for 
i — —N/2, . . . , N/2. The momenta {pi{0)} were drawn 
from a Maxwell-Boltzmann distribution at temperature 
T consistent with a given value of the energy density e. 
It is known that the dynamics of both the homogeneous 
(uniform mass, M ^ 1) FPU model [ill and the MFPU 
model 01 is strongly chaotic for large e values, whereas 
their dynamics, for small e values, corresponds to a chain 
of coupled harmonic oscillators. Consequently e was cho- 
sen in the range 0.01 < e < 100 with M = 1, 40, 60, 80, 
and 100. Finally, the 2N first-order Hamilton equations 
of motion were integrated using a third-order bilateral 
symplectic algorithm [l5| . which is a high-precision nu- 
merical scheme. 



III. STATISTICAL BEHAVIOR OF THE HEAVY 
IMPURITY 

Thermal equilibrium between the impurity and the 
FPU chain with N = 300 000 unit mass oscillators is 
attained within the time interval of t = 5 x 10^ nat- 
ural time units. Afterwards the heavy impurity per- 
forms Brownian motion for all e values studied |l^ . The 
momentum autocorrelation function (MACF) po{t) = 
{P{t)P{0))t/{P^iO))t of the heavy impurity was obtained 
by computing the time averages (• ■ •)( over a time inter- 
val oit = 2x 10^. For all e and M values considered the 
exponential fit exp(— </r), where r is the relaxation time 
of the MACF, is valid for t < 50. Since the magnitude of 
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FIG. 1: Relaxation time r vs energy density e. The continu- 
ous (dashed) lines correspond to the power-law fit e" for the 
e < 1 (e > 1) regime and all AI ^ 1 values. 



po(0 is negligible for t > 50 in all cases, its contribution 
was not considered in the computation of t. 

Fig-fflshows a graph of the relaxation time t vs the en- 
ergy density e in log-log scale for all the considered M val- 
ues, as e goes from the regular to the chaotic regime. We 
observe that the data points are separated into two dif- 
ferent and well-defined regions, depending on the e value. 
In all cases the dependence of r on the energy density is 
weak when e < 1. On the contrary, when e > 1, r de- 
creases rapidly as e increases. In each of these regimes 
T has a power-law scaling Tj^j{e) = ^,je°'M , with the 
same value for all the data points with the same M 
value. The slopes of each of the fits in Fig. ^ are approx- 
imately the same for the e < 1 regime. The same result 
occurs for e > 1, though with another slope value. These 
facts imply that there is a common Af-independent scal- 
ing exponent a for each regime. Following the method- 
ology of Ref. ^3 we estimate a ~ —0.019 for e < 1 and 
a K, —0.182 for e > 1. Thus we can conclude that the 
power-law scaling exponent undergoes a sudden change 
at a threshold value ec ~ 1.00, since the estimated a 
values differ by one order of magnitude. 



IV. PHASE-SPACE DYNAMICS 

In Ref. ^3 the change in the scaling behavior of r 
depicted in Fig. ^ was attributed to the change in the 
dynamical behavior of the chain as the energy density 
e goes from a regime of low chaos to a regime of fully 
developed chaos. However, it remains unexplained why 
the aforementioned change occurs at a precise value Cc of 
the energy density. To address this problem we consider 
in more detail the dynamical behavior of the FPU chain. 
As is known, the LLE Ai(e) is a parameter that can be 
used to quantify the degree of chaos in the dynamics of a 
given system with either low or large number of degrees 
of freedom. For the homogeneous FPU model [ij it is 
known that the behavior of the LLE for large (small) val- 
ues of the energy density e is Ai(e) t" . The exponent a 
has a large value in the weakly chaotic regime e < and 
a small value in the strongly chaotic regime e > e^, . That 
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FIG. 2: Largest Lyapunov exponent Ai vs energy density e 
for the homogeneous FPU model (M = 1). Asterisks corre- 
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FIG. 3: Largest Lyapunov exponent Ai vs energy density e for 
the MFPU model and all M values. References to power-law 
scalings valid for small and large e values are shown by dashed 
lines. Vertical dashed-dotted line indicates the approximate 
location of . Open symbols correspond to TV = 2000 and 
filled symbols to iV = 300 000. 



is, the exponent u undergoes a sudden change around the 
threshold e^. This provides an operational definition of 
the SST, which is defined by a crossover in the scaling 
behavior of Ai (e) . Further analytical studies [T^ . IT^ of 
the dependence of the LLE on the energy density e in 
homogeneous FPU lattices confirm the conjecture that 
the LLE reaches a finite, e-dependent, value in the ther- 
modynamic limit iV — > oo. 

For the homogeneous FPU model we have numeri- 
cally computed the LLE by a standard technique 
for chains with N = 300 000 and N = 2000 unit mass 
oscillators, as shown in Fig. |21 The results for these two 
N values overlap very well over the entire range of the e 
values studied with the best available anal ytic al estimate 
of Ai(e) valid in the thermodynamic limit '13l|. Since our 
numerical data are independent of the approximations 
made in Ref. to obtain the analytical estimate, they 
provide an independent corroboration of the stability of 
Ai(e) in the thermodynamic limit. This result justifies 
the study of the dynamics of the FPU chain with N val- 
ues smaller than those needed to compute the relaxation 
time T of the MACF of the heavy impurity. 

Fig. O shows the results for the MFPU model with 



FIG. 4: Spectrum of Lyapunov exponents for e = 26, A'^ = 
100, and all the considered M ^ 1 values. 



N = 2000 unit mass oscillators. For all the M values 
considered the scaling behavior of Ai remains unaltered 
by the presence of the heavy impurity embedded in the 
chain for the entire e value range considered. This re- 
sult indicates that the dynamics of the system is domi- 
nated by the FPU chain, with no contribution from the 
heavy impurity. A rough estimate of the crossover en- 
ergy density, which defines a common SST for both 
the homogeneous FPU and MFPU models, can be ob- 
tained from Fig.Oby the intersection of the straight lines 
that describe the low and high energy power-law asymp- 
totic behavior of Ai(e). The crossover between the two 
power laws ocurrs at « 1.18. A straightforward ob- 
servation is that the computed SST value, which is 
a distinctive characteristic of the microscopic dynamics 
of the entire system, is very close to the value €c ~ 1.00 
obtained from the crossover in the scaling behavior of the 
relaxation time t of the MACF, which is a property of 
the heavy impurity. This result suggests that the change 
in the scaling behavior of r exhibits macroscopically the 
existence of the SST . 

To characterize chaos in more detail we compute the 
Lyapunov spectrum (LS) of positive exponents {Ai} us- 
ing a standard method 0], as implemented in Ref. [l8l |. 
for chains of TV = 100 unit mass oscillators and several 
M values. The results are presented in Fig. 01 The lin- 
ear shape of the LS has been previously obtained for the 
homogeneous FPU model with e values much larger than 
the SST 0. We see that the LS is unaltered by the 
heavy impurity; that is, the LS is dominated by the posi- 
tive Lyapunov exponents corresponding to the FPU chain 
without the impurity. This type of behavior has been 
reported for the case of a heavy tracer particle in a two- 
dimensional molecular fluid, where the tracer and fluid 
particles are hard-disks undergoing elastic collisions 2(| , 
but not for the case of a continuous potential, like the 
one considered in this work. 

Another useful quantity to characterize chaos is the 
Kolmogorov-Sinai (KS) entropy /i^.^ which, for conserva- 
tive systems, can be written as /i^g = for all Ai > 0. 
The KS entropy describes the mean information produc- 
tion rate caused by all positive Lyapunov exponents along 
a trajectory in phase-space and therefore measures the 
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FIG. 5: KS entropy vs energy density for A'^ = 100. Vertical 
dashed-dotted line indicates the approximate crossover point 
between the regimes of weak and strong chaos, with dashed 
lines indicating the asymptotic power-laws valid for e -C 1 
and e » 1. 



degree of stochasticity. Results for a chain with = 100 
and two M values are reported in Fig.|Sl We observe that 
the KS entropy, just as the LLE, exhibits a change in its 
scaling behavior between two different and well-defined 
regimes, characterized by the scaling laws /i^,, ~ e^ °^ 
for e ^ 1 and h^^ ~ e"'^^ for e ^ 1. The crossover 
between both regimes occurs at an energy density value 
e that is closer to the threshold value Cc of the relaxation 
time T than the threshold of the LLE. Thus the KS 
entropy, just as the LLE, is a suitable quantity to probe 
the dynamics in the considered e value range. A similar 
behavior in the KS entropy has been previously reported 
for the one dimensional 0* model, with a threshold value 
of the energy that separates two integrable limits of this 
system plf. 



V. DISCUSSION 

From a macroscopic, i.e. statistical, perspective it is 
not at all clear what kind of collective mechanism, if 
any, could be held responsible for the crossover value 
Ec at which the scaling exponent a changes its value. 
However, from a microscopic Hamiltonian perspective 
the change of the scaling exponent a can be straight- 
forwardly interpreted as a macroscopic manifestation of 
the chaotic transition described by the SST. This inter- 
pretation arises by comparing Figs. ^ and |2| which sug- 
gest that the behavior of T„(e) and Ai(e) is due to the 
same mechanism. Since the microscopic dynamics of the 
MFPU and FPU models are practically identical, we can 
apply the same explanation of the origin of the chaotic 
transition that was previoulsy given in context of the 
FPU model The Hamiltonian of either the FPU or 
the MFPU model can be written in the form 

^eeM«i, (2) 

ll^oll 

where {9, 1) are the action-angle canonically conjugated 
variables and || • ■ • || is a suitable norm. A consequence 



of the perturbation is that the resonant manifolds 
n • = of iJ^ are destroyed for any small /x and are 
replaced by finite-thickness chaotic layers (n is an integer 
component vector and w is a vector whose components 
are LOi — dH^/dli). As this chaotic surfaces intersect 
the constant energy hypersurface for » 1, a chaotic 
network (the Arnold web) is produced which is every- 
where dense in phase space. At e > ec the resonances are 
strongly overlapped and microscopic diffusion is allowed 
in every direction in phase space, with a large value of 
the LLE. In this dynamical regime the scaling exponent 
a has also a large value, which produces a rapid decay 
in the value of the relaxation time r of the MACE of 
the heavy impurity. On the contrary, at e < ec, the res- 
onance overlapping is drastically reduced as the energy 
density decreases, diffusion in phase space occurs only 
along resonances, and the LLE takes a smaller value. 
Correspondingly, the scaling exponent a is smaller, which 
accounts for the almost constant value of the relaxation 
time r in this e value range. The behavior of the relax- 
ation time r(e) has a direct relationship with diffusion in 
configuration space, since r(e) is directly related to the 
self-diffusion coefficient of the heavy impurity through 
the Green-Kubo relation 

Although it is not our aim to directly address the prob- 
lem of the origin of diffusion in many-particle Hamilto- 
nian systems, a discussion will clarify some of our results. 
Recently a class of one-dimensional maps has been re- 
ported which present normal diffusive- like behavior in the 
absence of chaos (2^] . For the case of the MFPU model 
the diffusion coefficient of the heavy impurity obeys a 
unique power-law D ^ ^0.964 ^]^^ entire range of e val- 
ues considered 121. From these results it is clear that 
macroscopic diffusion is a combined result of the large 
number of degrees of freedom, the mass of the heavy 
impurity p^ . and the random initial conditions, rather 
than of the degree of chaos in the system. However, al- 
though diffusion itself is a collective phenomenon and, 
therefore, largely independent of the microscopic dynam- 
ics, some of its specific features, such as the crossover in 
the scaling behavior of r, can yield information on the 
dynamics of the system to which the heavy impurity is 
coupled, which in the case of the FPU chain studied here 
is the transition from weak- to a strongly chaotic regime. 
Nevertheless, in order to consider the heavy impurity as 
an effective probe of the Hamiltonian dynamics we have 
to study a suitable macroscopic variable, that in our case 
is the relaxation time r, which is computed from the mo- 
mentum or velocity of the heavy impurity. 

Finally we would want to remark some implications 
of the behavior of the KS entropy h^^. Although r is 
a property computed from the momentum of the heavy 
impurity, whereas the KS entropy is essentialy a property 
computed from the full dynamics of the FPU chain, both 
quantities provide information about the microscopic dy- 
namics. In fact, both reflect the chaotic transition at the 
SST e^, as shown in Figs. 13 and [S] The important point 
to be stressed is that the KS entropy is a property that 
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can be computed by other methods, such as time series- 
analysis [23. Recent theoretical work suggests that the 
information extracted from the time-series of the position 
of a Brownian particle is unable to unambiguously deter- 
mine the nature (either regular, chaotic, or stochastic) of 
the system to which the particle is coupled . This re- 
sult is entirely consistent with the lack of any signature of 
the chaotic transition associated to the SST in the behav- 
ior of the self-diffusion coefficient as a function of e in the 
case of our MFPU model pj]. In Ref. [H the KS entropy 
of an impurity in a harmonic crystal, which corresponds 
to our MFPU model for e ^ 1, was computed using the 
time-series of the position of the impurity. However, if 
the KS entropy is computed using the full dynamics of 
the system, as was done in the present paper for the case 
of the MFPU model, it can indeed detect the stochas- 
ticity transition, as is evident by inspecting Fig. [S] By 
comparing Figs. 13 and |S1 it can be infercd that Ai and 
hj^g convey almost the same information, which in turn 
is reflected in the behavior of the relaxation time T(e) 
as depicted in Fig. ^ Now, r(e) is computed from the 
MACF, which is a property that depends only on the 
momentum of the impurity. This observation suggests 
that, if the momentum time-series, instead of the posi- 
tion time-series, is used to compute the KS entropy by 
methods of time-series analysis, the stochasticity transi- 
tion associated with the SST of the FPU model could 
indeed be detected. We stress that we are not imply- 
ing that methods based on time-series analysis can de- 
tect chaos in a generic dynamical system. Even for one- 
dimensional systems, this possibility remains controver- 
tial |25|. The dynamics of three-dimensional systems or 
those with long-range interactions is still not well char- 
acterized, which poses an obstacle when trying to define 
which specific features of the microscopic dynamics of 
these type of systems could be detected at a macroscopic 



level. But, for the particular case of the MFPU model, 
it can be conjectured that the KS entropy, when com- 
puted from the momentum time-series, would display a 
different behavior at small e values compared to the cor- 
responding behavior at large e values. 



VI. CONCLUSIONS 

In this work we have presented evidence, obtained from 
a systematic study of certain dynamical parameters, that 
the inclusion of a heavy defect in a FPU chain does 
not affect its Hamiltonian dynamics. Furthermore, the 
crossover in the scaling behavior of both the LLE and the 
KS entropy, which depend on the dynamics of the whole 
system, has a macrospcopic manifestation in a similar be- 
havior of the relaxation time r of the MACF of the heavy 
impurity alone. An interesting development in this direc- 
tion would be to explore the possibility that certain dy- 
namical features could be detected in more complicated 
systems by monitoring the apropiate variables. Another 
open problem is to investigate if this transition between 
two diferent dynamical regimes can be detected by other 
methods, such as those provied by nonlinear time-series 
analysis. 
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